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Abstract. Extending the rigorous presentation of the "classical umbral calculus" 
|28) . the so-called partition polynomials are interpreted with the aim to point out 
the umbral nature of the Poisson random variables. Among the new umbrae intro- 
duced, the main tool is the partition umbra that leads also to a simple expression of 
the functional composition of the exponential power series. Moreover a new short 
proof of the Lagrange inversion formula is given. 



1 Introduction 

The symbolic method, nowadays known as umbral calculus, has been exten- 
sively used since the nineteenth century although the mathematical commu- 
nity was sceptic of it, maybe owing to its lack of foundation. This method 
was fully developed by Rev. John Blissard in a series of papers beginning 
from 1861 [HI-i-EH, nevertheless it is impossible to attribute the credit of the 
originary idea just to him since the Blissard's calculus has a mathematical 
source in the symbolic differentiation. In |22| Lucas even claimed that the 
umbral calculus has its historical roots in the writing of Leibniz for the suc- 
cessive derivatives of a product with two or several factors. Moreover Lucas 
held that this symbolic method had been subsequently developed by Laplace, 
by Vandermonde, by Herschel and augmented by the works of Cayley and 
of Sylvester in the theory of forms. Lucas's papers attracted considerable 
attention and the predominant contribution of Blissard to this method was 
kept in the background. Bell reviewed the whole subject in several papers, 
restoring the purport of the Blissard's idea 0] and in 1940 he tried to give 
a rigorous foundation of the mystery at the ground of the umbral calculus 
UJ. It was Gian-Carlo Rota who twenty-five years later disclosed the 
"umbral magic art" of lowering and raising exponents bringing to the light 
the underlying linear functional. In [221 and (21] the ideas from [25] led Rota 
and his collaborators to conceive a beautiful theory originating a large va- 
riety of applications. Some years later, Roman and Rota gave rigorous form 
to the umbral tricks in the setting of the Hopf algebra. On the other hand, 
as Rota himself has written |2B]: "...Although the notation of Hopf algebra 
satisfied the most ardent advocate of spic-and-span rigor, the translation of 
"classical" umbral calculus into the newly found rigorous language made the 
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method altogether unwieldy and unmanageable. Not only was the eerie feel- 
ing of witchcraft lost in the translation, but, after such a translation, the use 
of calculus to simplify computation and sharpen our intuition was lost by the 
wayside..." Thus in 1994 Rota and Taylor started a rigorous and simple 
presentation of the umbral calculus in the spirit of the founders. The present 
article refers to this last point of view. 

As it sometimes happens in the practice of the mathematical investigation, 
the subject we deal with does not develop the originary idea from which our 
research started in the spring of 1997, but this paper is closely related to 
it. In that period, Gian-Carlo Rota was visiting professor at the University 
of Basilicata and, during one of our latest conversations before his leaving, 
he shared with us his close interest for a research project: a combinatorial 
random variable theory. The delicate question arising from the underlying 
foundation side and the left short time led us to protract the discussion via e- 
mail intertwining it with different activities for several months. The following 
year, Gian-Carlo Rota held his last course in Cortona and we did not miss 
the opportunity to spend some time with him. We resumed the thread of our 
conversations and presented him with the doubts that gradually took hold 
of us. As usually, his contribution disclosed new horizons that have led us to 
write these pages. 

Our starting point is the umbral notion of the Bell numbers. Many classi- 
cal identities relating to these numbers are expressed in umbral notation at- 
taining up to a new umbra, the partition umbra, connected with the so-called 
"partition polynomials" generated by expanding the exponential function 
exp(/(x)) into an exponential power series. The hereafter developed theory 
of the Bell umbrae is not only an example of the computational power of 
the umbral calculus but it offers, we would like to believe, a natural way to 
interpret the functional composition of exponential power series tested by a 
new proof of the Lagrange inversion formula. Here the point operations ex- 
tended with a new one play a central role. From a probabilistic point of view, 
the functional composition of exponential power series is closely related to 
the family of Poisson random variables so that these random variables have 
found a natural umbral interpretation through the Bell umbrae. In partic- 
ular the probabilistic counterpoint of the partition umbra is the compound 
Poisson random variable. Also the less familiar randomized Poisson random 
variable unexpectedly find an umbral corresponding and in turn the umbral 
composition gives the way to generalize this last Poisson random variable. 

What's more, we believe that the probabilistic interpretation of the parti- 
tion umbra and of the umbral composition could give a probabilistic meaning 
to the Joyal species theory [2P| , namely a combinatorial random variable the- 
ory that we hope to deal in forthcoming publication. 
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2 The classical umbral calculus 

We take a step forward in the program of the rigorous foundation of the 
classical umbral calculus initiated by Rota and Taylor [57], |32|. 

In the following we denote by R a commutative integral domain whose 
quotient field is of characteristic zero and by A = {a, (3, . . .} a set whose 
elements are called umbrae. An umbral calculus is given when is assigned a 
linear functional E : R[A,x,y] — > R[x, y] such that: 



i) E[l] = 1; 

n) E[a i p ■■■ 1 k x n y rn ] = x n y m E[a i \E{^\ ■ ■ ■ E[-f k ] for any set of distinct 



umbrae in A and for k, n, m nonnegative integers (uncorrelation 

property); 

Hi) it exists an element e G A such that E[e n ] — 5o t n, for any nonnegative 
integer n, where 



iv) it exists an element u 6 A such that E[u n ] = 1, for any nonnegative 



The umbra e is named augmentation as Roman and Rota first called it |24| . 
We will call the umbra u the unity umbra. 

A sequence do, a%, a,2, ■ ■ ■ in R[x, y] is said to be umbrally represented by 
an umbra a when 



so that the linear functional E plays the role of an evaluation map. 

As Rota suggested, there is an analogy between umbrae and random vari- 
ables (r.v.) (see |32| ) . so we will refer to the elements in R[x, y] as moments 
of the umbra a. The umbra e can be view as the r.v. which takes the value 
with probability 1 and the umbra u as the r.v. which takes the value 1 with 
probability 1. 

An umbra is said to be a scalar umbra if the moments are elements of R 
while it is said to be a polynomial umbra if the moments are polynomials. 
Note that if the sequence do, <xi, ai, ■ ■ ■ is umbrally represented by a scalar 
umbra a, then it is ao = 1. In the same way, for polynomial umbrae, a 
sequence of polynomials po(x) , Pi(x) , P2(x) , . . . will always denote a sequence 
of polynomials with coefficients in R such that po(x) — 1 and p n [x) is of 
degree n for every positive integer n. 

A polynomial p £ R[A] is called an umbral polynomial. The support of p is 
defined to be the set of all occurring umbrae of A. Two umbral polynomials 
are said to be uncorrelated when their support are disjoint. 

If a and (3 are either scalar either polynomial umbrae, we will say that a 
and (3 are umbrally equivalent when 




integer n. 



E[a l ] = a 



for i = 0,1,2,... 



E[a]=E[0], 
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in symbols a ~ (3. Two scalar (or polynomial) umbrae are said to be similar 
when E[a k ] = E[/3 k ], k = 0, 1, 2, . . . or 

a fc ~/? fc , fc = 0,1,2,... 

in symbols a = /3. The notion of equivalence and similarity for umbral poly- 
nomials is obvious. 

The formal power series 

^ n\ 

n>l 

is said to be the generating function of the umbra a. Moreover, if the sequence 
ao, at, 0,2, • ■ ■ has (exponential) generating function f{t) and is umbrally rep- 
resented by an umbra a then i?[e Qt ] = /(<), in symbols e at ~ /(*)• When 
a is regarded as a r.v., fit) is the moment generating function. The notion 
of equivalence and similarity are extended coefficientwise to the generating 
functions of umbrae so that a = /3 if and only if e at ~ e' 3 *. Note that e e * ~ 1 



2.1 The point operations 

The notion of similarity among umbrae comes in handy in order to express 
sequences such 

53 ( { ) a i a n-n " = 0, 1, 2, . . . (1) 

i=0 ^ ' 

as moments of umbrae. The sequence Q cannot be represented by using only 
the umbra a with moments ao, at, a?, ■ ■ ■ because a{a n -i could not be written 
as E[a l a n ~ 1 ], being a related to itself. If we will assume that the umbral 
calculus we deal is saturated j2Z], the sequence ao,at,a2, ■ ■ ■ in R[x,y] is 
represented by infinitely many distinct (and thus similar) umbrae. Therefore, 
if we choose two similar umbrae a, a' , they are uncorrelated and 

= E[(a + a') n ]- 

Then the sequence QJ represents the moments of the umbra (a + a'). This 
matter was first explicitly pointed out by E.T. Bell (Sj who was not able to 
provide an effective notation: 

1 "...If in ax\ ■ ■ ■ -t-fiE there are precisely T summands ax, . . . , £x each of which 
is a scalar product of a scalar and x, we replace (— >) the T x's by T distinct umbrae, 
say a, ... ,x in any order, and indicate this replacement by writing 

ax+ ■ ■ ■ +£x — > aa+ ■ ■ ■ +£x. 
1 The quotation needs more details. It is 

ax\ Vex = (ax H h ex , ax N H h ex N , . . .), 
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Then (aa+ ■ ■ ■ +£x) is to be calculated by (1.22) and the exponents are degraded 
as in (1.20). In the result, each of a, . . . , x is replaced (<—) by x; the resulting poly- 
nomial is defined to be the N— th power (otx+ ■ • • +£x) N of the sum (ax+ ■ • ■ +£x). 
For example 

(ax + /3x) 3 — > (aa + /3x) 3 ; 

(aa + f3x) :i = a 3 a3Xo + 3a 2 /3a,2Xi + 3a/3 2 aiX2 + (3 3 aoX3, 
<— a 3 X3Xo + 3a 2 (3x2Xi + 3af3 2 x\X2 + (3 3 xoX3; 
(ax + (3x) 3 = (a 3 + /3 3 )x x 3 + 3af3(a + /3)x 1 X2..." 



The last identity makes sense when the left side is replaced by (ax + fix 1 ) 3 
with x = x' , but Bell did not have the notion of similar umbrae. However, the 
need of handling sequences like Q leads to introduce some new operations 
between umbrae, as showed in the next sections. 



2.2 The point product 

We shall denote by the symbol n.a an auxiliary umbra similar to the sum 
a' + a" + . . . + a'" where a', a", ... , a'" are a set of n distinct umbrae each 
of which is similar to the umbra a. We assume that O.a is an umbra similar 
to the augmentation e. A similar notion n.p is introduced for any umbral 
polynomial. The following statements are easily to be proved: 

Proposition 1. (i) If n.a = n./3 for some integer n ^ then a = /3; 

(ii) if c G R then n.(ca) = c(n.a) for any nonnegative integer n; 

(Hi) n.(m.a) = (nm).a = m.(n.a) for any two nonnegative integers n, to; 

(iv) (n + m). a = n.a + to. a' for any two nonnegative integers n, m and any 
two distinct umbrae a = a'; 

(v) (n.a+n.p) = n.(a+/3) for any nonnegative integer n and any two distinct 
umbrae a and (3. 

Proposition 2. If a is an umbra with generating function e at ~ f(t), then 
the umbra n.a has generating function e' n Q )* ~ [f(t)\ n . 

Proof. It follows from the definition of the auxiliary umbra n.a. □ 
where x is an umbra. The formula (1.22) is 

(aa+ ■ ■ ■ +£x) N = V M Sl ,...,s T a Sl ■ ■ ■ £ ST a Sl ■ ■ ■ x St 

with Msi^.^Sr the coefficient of xf 1 ■ ■ ■ x^ T in the expansion of (x\ + ■ ■ ■ +xt) n 
through the multinomial theorem. The formula (1.20) is 

(aa-\ Vex) N = ^Ms lr ..,s T Q Sl • • ■ e St as x ■ ■ -xs T . 
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The moments of the umbra n.a are the following polynomials in the variable 
n 

k 

E[(n.a) k }=q k (n) = J2( n ^ B ^ £ = 0,1,2,... (2) 

i=0 

where Bk t i = -Bfe,i(ai, a2, ■ ■ ■ , Ofc-i+i) for i < k are the (partial) Bell expo- 
nential polynomials 2 , (n)i is the lower factorial and are the moments of 
the umbra a. Recalling that 

oo 1 

=*[/(*)-!]'. ( 3 ) 



the identity J2J follows from 



[/(or - Ew.M = E ( £(»)<**.<] t (4) 



fc=0 \i=0 



If in (|2J| set a = it, then gfc(n) = n fe . Note that go ( n ) = 1> ?fc(0) = and the 
polynomial sequence {qk(n)} is of binomial type as it follows by using the 
statement (iv) of Proposition 1: 

[(n + m).a] k ~ [n.a + m.a'] k — £ ( • ) ( m - a ') fc 

»=o ^ ' 

Moreover a variety of combinatorial identities could be umbrally interpreted. 
As instance in point, the classical Abel identity becomes 

(a + (3) n ^J2(l) a ( a - k ^) k ' 1 ((3 + k.y) n - k , n= 0,1,2,.... (5) 



k=Q 



The expression of the polynomial sequence {qk(n)} in (J2J suggests a way 
to define the auxiliary umbra x.a when x £ R, however it is impossibile to 
give an intrinsic definition. Up to similarity, the umbra x.a is the polynomial 
umbra with moments 

k 

E[{x.a) k ]=q k {x)=Y,^ Bk ^ k = M,2,.... (6) 



i=0 



Note that qk(x) = x k when a = u. 

Proposition 3. If a is an umbra with generating function e at ~ f(t), then 
the umbra x.a has generating function e^")' ~ [f(t)] x . 

Proof. It follows from (JIJ and J21 with n replaced by x. □ 

Corollary 1. (i) If x.a = x.j3 for x 6 R — {0} then a = (3; 
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(ii) if c 6 R then x.(ca) = c(x.a) for any x £ R; 
(Hi) x.{y.a) = {xy).a = y.(x.a) for any x,y £ R; 

(iv) (x + y).a = x.a + y.a' for any x,y £ R and any two distinct umbrae 
a = a'; 

(v) (x.a + x.p) = x.(a + (3) for any x € R and any two distinct umbrae a 
and (3. 

Theorem 1. Up to similarity, each polynomial sequence of binomial type is 
umbrally represented by an auxiliary umbra x.a and viceversa. 

Proof. ^From the statement (iv) of the corollary 1, it follows that the polyno- 
mial sequence {qk{x)} is of binomial type. Viceversa, first observe that from 
® it is 

D x [q k {,x)] x= n = a k + F(a 1 ,a 2 , . . .,a k -i) (7) 

where F is a function of the moments a±, a 2 , . . . , ak-i- Let {pk(x)} be a poly- 
nomial sequence of binomial type. Through (JJJ, the moments of the umbra a 
are uniquely determined by the knowledge of the first derivative respect to x 
of pk(x) evaluated in 0. Moreover, the sequences of first derivate respect to x 
of Pk(x) evaluated in uniquely determines a sequence of binomial type. □ 

Similarly with it has been done for the auxiliary umbra x.a, we define a 
point product among umbrae. Up to similarity, the umbra (3.a is an auxiliary 
umbra whose moments are umbrally expressed through the umbral polyno- 
mials q a ,k(P) ■ 

fe 

{iS.af ~g Q , fc (/3) =J2(@)iBk,i k = 0,1,2,.... (8) 

If a is an umbra with generating function e Qt ~ f(t), then the identity 
could be rewritten as 

i=0 k=0 \i—0 / 

so that e^' a '* ~ [f(t)]&. Moreover if (3 is an umbra with generating function 
e' 3 * ~g(t), then 

[/( t )]/3~ e W(t)~ ff [log/(t)]. 
This proves the following proposition. 

Proposition 4. If a is an umbra with generating function e at ~ f(t) and 
(3 is an umbra with generating function e^ 1 ~ g(t), then the umbra (3. a has 
generating function 

eC-)* ci[f(t)f~ g [log f(t)}. (10) 
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Corollary 2. If 7 = 7' then 

(a + (3)."f = a.-f + (3.J 1 . 
Proof. Let e 7 * ~ h(t) the generating function of the umbra 7. It is 

e [(a+/J).7]t ^ [^(t)]^/ 3 ~ [h{t)] a [h(t)f ~ e^^^e^' 7 ')* 
from which the result follows. □ 

Remark 1. As Taylor suggests in |32j . the auxiliary umbra /3.a provides an 
umbral interpretation of the random sum since the moment generating func- 
tion g[\og f(t)} corresponds to the r.v. Sn = X\ + X 2 + ■ ■ ■ + Xn with Xi in- 
dipendent identically distributed (i.i.d.) r.v. having moment generating func- 
tion f(t) and with N a discrete r.v. having moment generating function g(t). 
The probabilistic interpretation of the corollary 2 states that the random sum 
Sn+m is similar to S n + Sm, where N and M are two indipendent discrete 
r.v. 

The left distributive property of the point product respect to the sum 
does not hold since 

e M/3+7)]* ~ [g(tW[h(t)F gk f[logg(t)]f[logh(t)] 

where g(t) ~ e^, h(t) ~ e 7 * and /(*) ~ e at . 

Again this result runs in parallel with the probability theory. In fact, let 
Z = X + Y be a r.v. with X and Y two indipendent r.v. The random sum 
Sn — Z\ + Zi + • ■ • + Zn, with Zi i.i.d. r.v. similar to Z, is not similar to the 
r.v. + Sjf where — X\ + X 2 + • • • + Xjy and Xi i.i.d. r.v. similar to 
X, and where S 1 ^ = Y\ + Y% + ■ ■ ■ + Y/v and ^ i.i.d. r.v. similar to Y. 

Corollary 3. (i) If [3. a = (3.-f then a = 7; 

(ii) if c 6 R then (3.(ca) = c(/3.a) for any two distinct umbrae a and (3] 
(in) (3.{^.a) = (p.j).a. 

Proof. Via generating functions. □ 

To end this section, we deal with the notion of the inverse of an umbra. Two 
umbrae a and (3 are said to be inverse to each other when a + (3 = e. Recall 
that, dealing with a saturated umbral calculus, the inverse of an umbra is 
not unique, but any two inverse umbrae of the umbra a are similar. 

Proposition 5. If a is an umbra with generating function e at ~ f(t) then 
its inverse (3 has generating function e^' ~ 

Proof. The result follows observing that e' Q+ ' 3 - ) * ~ 1. □ 
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Similarly, for every positive integer n and for every umbra a € A, the inverse 
of the auxiliary umbra n.a, written as — n.a' with a = a', is similar to 

f3' + f3" H h [3'" where f3' , f3" , . . . , (3'" is any set of n distinct umbrae similar 

to (3, being (3 the inverse of a. The notation —n.a' is justified by noting that 

n.a — n.a' = (n — n).a = O.a = e. 



Proposition 6. If a is an umbra with generating function e at ~ f(t), then 
the inverse of n.a has generating function e { -~ n - a ~ [f(t)]~ n . 

Proof. The result follows observing that e ( n - a ~ n - a )* ~ 1. □ 

The inverse of the umbra x.a is the umbra —x.a' with a = a' defined by 

x.a — x.a = (x — x).a = O.a = e. 

2.3 The point power 

As it is easy to be expected, the definition of the power of moments requires 
the use of similar umbrae and so of a point operation. This notion comes 
into this picture by a natural way, providing also an useful tool for umbral 
manipulation of generating function. 

We shall denote by the symbol or™ an auxiliary umbra similar to the 
product a' a" ■ ■ ■ a'" where a', a" , . . . , a'" are a set of n distinct umbrae each 
of which is similar to the umbra a. We assume that a is an umbra similar to 
the unity umbra u. A similar notion is introduced for any umbral polynomial 
p. The following statements are easily to be proved: 

Proposition 7. (i) If c e R then (ca)- n = c"ar" for any nonnegative integer 
n^O; 

(ii) (a- n )- m = a- nm = (cr m )'" for any two nonnegative integers n,m; 

(Hi) ar(" +m ) = a- n (a')- m for any two nonnegative integers n, m and any two 

distinct umbrae a = a'; 
iv) (ar™) fe = (a )"" for any two nonnegative integers n, k. 

By the last statement, the moments of or™ for any integer n are: 

E[{a- n ) k ] - E[{a k )- n ] = a n k , k = 0,1,2,... (11) 

so that the moments of the umbra or™ are the n— th power of the moments 
of the umbra a. 

Proposition 8. The generating function of the n—th point power of the um- 
bra a is the n—th power of the generating function of the umbra a. 



10 E. Di Nardo, D. Senato 



Note that, by virtue of Propositions 2 and 8 it is 

e (n. a )t _ ( e at).n ( 12 ) 

The relation (|12|) restores the natural umbral interpretation of [f(t)] n . More, 
let us observe that if a and /3 are not similar, it is 



,io.(n-i) 



The point power operation leads us to define the point exponential of an 
umbra. We shall denote by the symbol e. Q the auxiliary umbra 



a 

e." =u ' 

n 

n=l 



We have immediately e. e = u. 
Proposition 9. For any umbra a, it is 

e .(n.a) ^ ( e «)•*». (14) 



Proof. It results 

£[(e a )-"] = £[e Q ]" = e nE[a] = ^ 

and also 



n k E\oe] 



k\ 

k>0 



mp („.«)! _ \- £[(rc.a)- fc ] _ \- n k E[a k ] 
1 k\ ~ ^ k\ ' 

fc>0 fe>0 

by which (fTljl follows. □ 

Up to similarity, the expression of the moments given in l|ll|l justifies the 
definition of the auxiliary umbra a x as the umbra whose moments are 

E[(a-*) k ] = a%, k = 0,1,2,.... 



Proposition 10. Let a be an umbra and e at ~ f(t) its generating function. 
It is 

e (x.a)t _ ^atyx _ [f(f^. 

Via moments, it is possible to prove the analogue of Proposition 7 where n 
and m are replaced by x and y with x,y £ R. 

Once again, we define the auxiliary umbra cc? as the umbra whose mo- 
ments are umbrally equivalent to 

(or^-af, k= 0,1,2,... 

and we set e a = e. 
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Proposition 11. (i) (or' 3 )' 7 = a-' 7 '' 3 '; 

(ii) a'' /3+7 ' = ofP(a')^ for any two distinct umbrae a = a'. 
Proof. It follows via moments. □ 

Proposition 12. Let e at ~ f(t) be the generating function of the umbra a. 
It is 

In closing, we notice that the generating function of the point product be- 
tween umbrae is umbrally equivalent to the following series: 

00 \p at 

eW'^^li-l (15) 

i=0 

by the relation and Proposition 8. 
3 Bell umbrae 

The Bell numbers B n have a long history and their origin is unknown: Bell 
ascribes them to Euler even without a specific reference Usually they 
are referred as the number of the partitions of a finite nonempty set with n 
elements or as the coefficients of the Taylor series expansion of the function 
exp(e* — 1). It is just writing about the Bell numbers that Gian-Carlo Rota 
gives the first glimmering of the effectiveness of the umbral calculus in 
manipulating number sequences, indeed his proof of the Dobinski's formula 
is implicitly of umbral nature. 

In this section, the umbral definition of the Bell numbers allows the proofs 
of several classical identities (cf. j^H]) through elementary arguments and 
smooths the way to the umbral interpretation of the Poisson random vari- 
ables. 

Definition 1. An umbra j3 is said to be a Bell scalar umbra if 

(/3)„~1 n = 0,l,2,... 

where (f3)o = 1 and (/?)„ = [3(f3 — 1) •••(/? — n + 1) is the lower factorial. 

Up to similarity, the Bell number sequence is umbrally represented by the 
Bell scalar umbra. Indeed, being 

n 

where S(n, k) are the Stirling numbers of second kind, then 



n n 

E((3 n ) = ]T S(n, k)E[(j3) k ] = S(n, k) = B, 

k=0 k=Q 



12 E. Di Nardo, D. Senato 

where B n are the Bell numbers. 

The following theorem provides a characterization of the Bell umbra. 

Theorem 2. A scalar umbra (3 is a Bell umbra iff 

j3 n+1 ~ ((3 + u) n n = 0,1,2,... (16) 

Proof. If (3 is the Bell scalar umbra, being (3((3 — u) n ~ (f3) n+ i it is E[(3{(3 — 
u) n ] = 1 = E[((3) n ]. By the linearity it follows 

E[(3p(f3 - «)] = S[p(/3)] 

for every polynomial p in j3. So, the idenity Ijl6|l follows setting p{(3) — ((3 + 
u) n . Viceversa, the relation l|16|l gives 

E[(3 n+1 ] =E[(f3 + ur\ =J2( n k )E[P k } 



fe=0 



or setting E[(3 n ] — B n one has 



- E u ) B 



k 

that is the recursion formula of the Bell numbers. □ 

Corollary 4. // (3 is the Bell scalar umbra, then 

D t [e^} ~ e W. (17) 

Proposition 13. // /3 is the Bell scalar umbra, for any integer k > and 
for any polynomial p(x) the following relation holds 

P{f3 + k.u)~{f3) k p{f3)~p{[3). 

Proof. For n > fc, by the definition 1 it follows 

{(3) n ~{(3) n+k ~{l3) k {(3-k.u) n . 

Thus for any polynomial q it is 

q{(3)~((3) k q{(3-k.u) 

by which one has 

(f3 + k.u) n ~(f3) k (3 n , n = 0,1,2,... 
setting q(/3) = (/3 + k.u) n . The result follows by linearity. □ 
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Proposition 14. The generating function of the Bell umbra is 



—u 



(18) 



Proof. By the definition 1 and the relation (|15|l it is 




Thus (TSJ follows from the relation JT3|). □ 

Remark 2. Let us go on with our probabilistic counterpoint noting that the 
Bell umbra can be view as a Poisson r.v. with parameter A = 1. Indeed, the 
moment generating function of the Bell umbra is exp(e* — 1) (see (^SJ) so 
that P(e*) = exp(e* — 1) where P(t) is the probability generating function, 
and therefore P(s) = exp(s — 1). By that, the moments of a Poisson r.v. with 
parameter 1 are the Bell numbers and its factorial moments are equal to 1. 

The following theorem makes clear how the proof of Dobinski's formula be- 
comes natural through the umbral expression of generating function. 

Theorem 3 (Umbral Dobinski's formula). The Bell umbra (5 satisfies 
the following formula: 



Proof. Being e^* ~ e. u e. e "* it is 

00 Jk.u)t °° ~\ ( 00 (h n\ n t n 1 

fc=0 k=Q (n=0 J 

by which the result follows. □ 

3.1 The Bell polynomial umbra 

Definition 2. An umbra <f> is said to be a Bell polynomial umbra if 




(</>)„ ~ x 



n = 0,l,2,.... 



Note that <f> = for x = 1. Moreover, being 



k=0 



by the definition 2 it follows 



7 ; 



n 



(19) 
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The polynomials <P n (x) have a statistical origin and are known in the liter- 
ature as exponential polynomials. Indeed, they were first introduced by Stef- 
fensen [2] and studied further by Touchard [233 an d others. Rota, Kahaner 
and Odlyzko |26| state their basic properties via umbral operators. 

Proposition 15. The generating function of the Bell polynomial umbra is 



Thus J20J) follows from the relation JT3|). □ 

The following theorem provides a characterization of the Bell polynomial 
umbra. 

Theorem 4. An umbra is the Bell polynomial umbra iff 



where [3 is the Bell scalar umbra. 
Proof. The result comes via l)20[l. □ 

Remark 3. The Bell polynomial umbra can be view as a Poisson r.v. with 
parameter A — x. Indeed, the moment generating function of the Bell poly- 
nomial umbra is exp[a;(e* — 1)] (see (jSOl) so that P{e l ) = cxp[a;(e i — 1)], where 
P(t) is the probability generating function, and therefore P(s) — exp[x(s— 1)]. 
By that, the moments of a Poisson r.v. with parameter x are the exponential 
polynomials and its factorial moments are equal to x n . 

When x — n, the Bell polynomial umbra n.(3 is the sum of n similar uncor- 
rected Bell scalar umbrae, likewise in probability theory where a Poisson r.v. 
of parameter n can be view as the sum of n i.i.d. (eventually uncorrelated) 
Poisson r.v. with parameter 1. More in general, the closure under convolu- 
tion of the Poisson probability distributions i.e. F s * F t — F s+t , where F t is 
a Poisson probability distribution depending on the parameter t, is umbrally 
translated by x.j3 + y.f3' = (x + y).(3 (cf. statement (iv) of Proposition 1). 

The next theorem is the polynomial analogue of the theorem. 
Theorem 5. An auxiliary umbra x.(3 is a Bell polynomial umbra iff 




(20) 



Proof. By the definition 2 and the relation (|15fl it is 




(f> — x./3 



(x.p) n+1 ~x(x.[3 + u) n , n = 0,l,2,. 



(21) 
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Proof. Observe that D^e^*] ~ D t [(e^)- X ] ~ xe^^^D^' 1 }, where 
0' = 13. From JTU) it is D^e^*] ~ xe^^^* from which the result follows 
immediately. Viceversa, the relation (|21(1 gives (|16|l for x = 1, by which it 
follows that /? is the Bell scalar umbra. □ 

The formula H21(l represents the umbral equivalent of the well known recursive 
formula for the exponential polynomials: 

k=0 ^ ' 

Similarly, the next proposition gives an umbral analogue of the Rodrigues 
formula for the exponential polynomials (cf. |2fi|V 

Proposition 16. The Bell polynomial umbra x.(3 has the following property: 

D x [(x.0) n }~(x.p + u) n -(x.p) n . 
Proof. ^From fl"%|) it is 

by which the result follows immediately. □ 

In closing we state the polynomial version of the umbral Dobinski's formula. 



Proposition 17. The Bell polynomial umbra x.j3 satisfies the following re- 
lation: 

OO / -J \fT fc 

(x.f3) n ~e-*- u E fc! X ■ 

fc=0 



Proof. Being e^* ~ er x - u eS x - eU ^ it i 



is 



llt\.fc ryk (k.u)t 

fc=0 ' fe=0 



fc! ^ fc! 

) 

by which the result follows. □ 



3.2 The exponential umbral polynomials 

Let us introduce a new family of umbral polynomials that turns out to be an 
useful tool in the umbral composition, also disclosing an unexpected proba- 
bilistic interpretation. 
Set 

n 

<P n (a)=Y,S{n,k)a k , n = 0,1,2,..., (22) 

fc=0 
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we will call <? n (a) exponential umbral polynomials. By the identity 1)19(1 . being 
(x.f3) n ~ # n (x), where (3 is the Bell scalar umbra, it is glaring that 

- («-/9) n n = 0,1,2,... 

and 

(a./?) n ~a n n = 0,1,2,..., (23) 

a formal proof passing through similar arguments already produced for the 
umbra x.(3. 

Proposition 18. Let (3 be the Bell scalar umbra. If e at ~ f(t) is the gener- 
ating function of the umbra a then 

e (a./3)* ~ yy _ !]. (24) 

Proof. The result follows by the relation I|1U|I observing that e^* ~ e e * _1 . □ 

When /(t) is considered as the moment generating function of a r.v. X, a 
probabilistic interpretation of (|2~l}l suggests that the umbra a. (3 represents a 
Poisson r.v. Nx with random parameter X. Indeed the probability generating 
function of Nx is 

P(s) = J2 P ( N x = k)s k = / P(N X = k\X = x)dF x {x) 

= E h / xke ~ XdF *w = E - l ^ Le ^ = ^ s x ) 

fe=0 ' ^° fc=0 

hence the moment generating function of A/jf is f(e t — 1). To the best of 
our knowledge, this r.v. has been introduced in [20] as the randomized Pois- 
son r.v. Once more the closure under convolution of the Poisson probability 
distributions leads us to claim that the point product a. (3 is the umbral cor- 
responding of the random sum of independent Poisson r.v. with parameter 1 
indexed by an integer r.v. X. 

4 The partition umbra 

As suggested in [2EJ , there is a connection between polynomials of binomial 
type and compound Poisson processes. Two different approaches can be found 
in |3U] and in [T2|. In this section, we suggest a way, that we believe to be 
natural, in order to make clear this connection. 

Definition 3. An umbra ip is said to be an a—partition umbra if 

ip = (3.a 

with (3 the Bell scalar umbra. 
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Note that the u— partition umbra is the Bell scalar umbra. 
Proposition 19. The generating function of the a— partition umbra ip is 

e^ ~ e. (eQt - u) . (25) 
Proof. ^From (|15() and by the definition 1 one has 



n=0 



E 



Thus follows from the relation 1(15)1 . □ 

The generating function (|25[1 leads us to interpret a partition umbra 
as a compound Poisson r.v. with parameter 1. As well known (cf. |19p. a 
compound Poisson r.v. with parameter x is introduced as a random sum 
Sm = X\ + X2 + ■ ■ ■ + Xm where N has a Poisson distribution with pa- 
rameter x. The point product of definition 3 fits perfectly this probabilistic 
notion taking into consideration that the Bell scalar umbra (3 plays the role 
of a Poisson r.v. with parameter 1. What's more, since the Poisson r.v. with 
parameter x is umbrally represented by the Bell polynomial umbra x./3, & 
compound Poisson r.v. with parameter x is represented by the polynomial 
a— partition umbra x.tp = x./3.a with generating function 

e (*-V)* ~ e [x.(e at -u)] t (26) 

The name "partition umbra" has also a probabilistic ground. Indeed the 
parameter of a Poisson r.v. is usually denoted by x = At, with t representing 
a time interval, so that when this interval is partitioned into non-overlapping 
ones, their contributions are stochastic independent and add to Sn- The last 
circumstance is umbrally expressed by the relation 

(x + y).f3.a = x./3.a + y.f3.a (27) 

that also assures the binomial property for the polynomial sequence defined 
by x.fi.a. In terms of generating functions, the formula (|27|l means that 

h x+y {t) = h x {t)h y (t) (28) 

where h x {t) is the generating function of x.(3.a. Viceversa every generating 
function h x {t) satisfying the equality is the generating function of a 
polynomial a— partition umbra, namely h x (t) has an umbral expression of 
the form i|2l)jl . 

Going back the moments of a partition umbra, according to the definition 
of the Bell scalar umbra and from |(HJ| it is 

n 

E[{(3.a) n ] = ^B nyk (a 1 ,a 2 , ■ ■ ■ ,a n - k+1 ) = F n (ai,a 2 , . . . ,a„) (29) 
fc=i 
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where Y n = Y n (a\, a 2 , . . . , a n ) are the partition polynomials (or complete Bell 
exponential polynomials) and Oi are the moments of the umbra a. Although 
the complexity of the partition polynomial expression, their umbral inter- 
pretation l|29[) allows an easy proof that they are of binomial type, simply 
observing that (3. a + /3.7 = {3. (a + 7). 

Partition polynomials have been first introduced by Bell PP who gave a 
pioneer umbral version of them in 2 . Because of their generality, they include 
a variety of other polynomials such as the cycle indicator of the symmetric 
group and other of interest in number theory. 

As already done for the Bell scalar umbra, the next theorem characterizes 
the partition umbrae and also provides the following recursive formula for the 
partition polynomials: 



^+1(01,02, • • • , o n +i) = y ][ . ) a„_fc+iYfc(ai,02, • • • , a k ). 



Theorem 6. Every a— partition umbra verifies the following relation 

{/3.a) n+1 ~ a'(f3.a + a') n a' = a, n = 0, 1, 2, . . . (30) 
and viceversa. 

Proof. Let tft an a— partition umbra. Then from l|25|) it is Dtle^'"^ 1 } ~ 
e (/3 - Q)t A[e Q '*], where a' = a. The identity (j20Jl follows observing that Dt [e Q t ] : 
a'e a *. Going back the previous steps, from (|30(l one has that (3. a has gener- 
ating function (|25(l and so it is an a— partition umbra. □ 

The moments of the polynomial a— partition umbra are 

n 

E[(x.f3.a) n ] = ^{x.(3) k B nyk (a 1 ,a 2 , a n - k +i) 

k=l 
a 

= ^ x k Bn tk (ai,a 2 , • • ■ , a„_ fc+ i) (31) 



fc=i 



according to the definition 2. The same arguments given in the proof of the 
theorem 6 lead to state that every polynomial partition umbra verifies the 
following formula 

(x.(3.a) n+1 ~ xa\x.(3.a + a') n a' = a, n = 0, 1, 2, . . . 



and viceversa. 
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4.1 Umbral expression of the functional composition 

An umbral wording of the functional composition of exponential formal power 
series is a thorny matter. It was broached by Rota, Shen and Taylor in 28 
passing through the sequence of Abel polynomials. In this last section, we give 
an intrinsic umbral expression of this operation via the notion of partition 
umbra. 

Definition 4. A composition umbra of the umbrae a and 7 is the umbra 

X = l-P-at 

where (3 is the Bell scalar umbra. 

In other words, the composition umbra \ is the point product of the umbra 
7 and the a— partition umbra [3. a. 

Remark 4- As already stressed in section 3.2, the umbra j.(3 represents a ran- 
domized Poisson r.v.. Hence it is natural to look at the composition umbra as 
a new r.v. that we will call compound randomized Poisson r.v. Moreover, be- 
ing (7./3).a = 7.(/3.a) (cf. statement (ii) of corollary 3), the previous relation 
allows to see this new r.v. from another side: the umbra 7.(/3.a) generalizes 
the concept of a random sum of i.i.d. compound Poisson r.v. with parameter 
1 indexed by an integer r.v. X, i.e. a randomized compound Poisson r.v. with 
random parameter X. 

Proposition 20. The generating function of the composition umbra ^f. (3. a is 
the functional composition of the generating functions e at ~ f(t) and e 7 * ~ 
<?(*)• 

Proof. Via (23 it is e ^ Q )* ~ e f{ ^~ x . The result follows by (HUJ observing 
that e W- a >]* ~ g {logle'W- 1 ]} . □ 

The moments of the composition umbra are 

n 

(a\, a 2 , ■ ■ ■ , a n -k+i) (32) 

k=a 

where aj are the moments of the umbra a. Indeed, by (JSJ it is 

n 

(~f.f3.aT ^ ]T(7-/?) fc Bn, fe 

fc=0 

and (O follows from (t53|) . 

Once more, we give a characterization of the composition umbra in the 
next theorem. 
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Theorem 7. Every composition umbra verifies the following relation 

{ 1 .(3.a) n+1 = ya'i'y.p.a + a') n a = oi, n = 0, 1, 2, . . . (33) 
and viceversa. 

Proof. Let \ a composition umbra of a and 7. Then from Proposition 20, 
it is A[e**] ~ g'[f(t) - l]f(t). Equation (gSJ) follows being f'(t) ~ a'e a>t 
with a' = a and g'[f(t) — 1] ~ je xt . Going back the previous steps, from 
l|33|) it follows that 7./?. a has generating function g[f(t) — 1] and so it is a 
composition umbra of a and 7. □ 

At this point, as custom, we put to test the definition 4 of composition umbra, 
giving a proof of Lagrange inversion formula. In the literature (cf. 18 for a 
plenty of references) different forms of the Lagrange inversion formula are 
derived using umbral calculus. The main tool of our proof is the umbral 
expression of the (partial) Bell exponential polynomials that we state in the 
next proposition. 

Lemma 1. It is 

B 1hk ( ai ,a 2 ,...,a n - k+1 )~ f"J a- k (k.a) n - k (34) 

where a is the umbra with moments E [a"l = — " +1 — r, n = 1, 2, . . . . 

ai(n + 1) 

Proof. By the identity J2J it results 

Bn,fc(oi,o 2s ... ) a B _fc + i) = ^D { t n) [{f(t) - l) fc ] t=0 

where D\ [-]{ = q is the n— th derivative respect to t evaluated in t = and 
f{t) ~ e at , so that 

B n .k{ai,a 2 , ■ ■ ■ ,a n _fc+i) ~ T,D^[(e at - w)' fc ] t= o- 

fc! 

On the other hand, by the moment expression of umbra a it follows e at — u ~ 
ate at . Therefore one has 

D^[(e at - u)- k ] ~ a- k D { t n) [t k e k ^} 

k 



a 

3=0 



j)Di f) [t k ]D^[e k -% 



using the binomial property of the derivative operator. Finally, the result 
follows evaluating the right hand side of the previous formula in t = and 
observing that D { t n ~ k) [e fc -"] t=0 - (k.a) n ~ k . □ 
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Remark 5. Let a = u. Then or fc ~ 1 for k = 0, 1, 2, . . . and i? n ^(1, 1, . . . , 1) = 
S(n,k) the Stirling number of the second kind. Moreover it is a = (—1.(5) 
where S is the Bernoulli umbra whose moments are the Bernoulli numbers 
(cf. [2HI)- From Lemma 1 it results 

S(n,k)^( n k )(-k.Sr- k 

as already stated by Rota and Taylor through a different approach (cf. Propo- 
sition 9.1 E7]) 

Theorem 8 (Lagrange inversion formula). Let e at ~ f(t) and e 7 * ~ 
<?(*)■ Ifg[f(t) - 1] - f[9(t) - 1] = 1 + 1 then 

a- k j k ~ (-k.a)^ 1 ^^ 1,2,.... (35) 
Proo/. By the formulas (jS2J and lj3H) , it is 

^^E(fc) a ' fc7fc( " fc - S)n_fe - (36) 

On the other hand, the Abel identity (JSJ) gives 

X" ~ E (fc) X(X - fc.a) fc - 1 (fc.a) n - fe . (37) 
fc=0 ^ ' 

Comparing (|36() with i|37|) one has 

by which the result follows expanding the right hand side of the previous 
formula by the binomial theorem and observing that from g[f(t) — 1] = 1 + t 
itisx-land X J -0,j = 2,3,.... □ 

More explicitly, the formula (I35|) says that the k— th coefficient of the gen- 
erating function g(t) is equal to the (k — 1)— th coefficient of the generating 
function [(/(£) -l)/i]~ fc , when £/[/(*) -1] = 1+t. Note that if /(t) — 1 = te~* 
then ar fe ~ 1, a = —l.u and from P5|l it is j k ~ (k.u) k ^ 1 . 

In closing, let us observe that if a\ = 1 then /(f) — 1 ~ te at and the 
Lagrange inversion formula (|35|l becomes 

7 fe = (-k.a) k -\ 

On the other hand, if the generating function g(t) is written as g(t) — 1 ~ ie 7 ' 
then the Lagrange inversion formula (|35|l becomes 

Iff-l ~ (-fc.S)^ 1 

that is equivalent to the version given in [2Hj by using the Abel polynomial 
sequence and its delta operator. 
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